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PROPAGATION  OF  SHORT  LASER  PULSES  IN  PLASMA  CHANNELS 


1.  Introduction 

Advances  in  laser  technology  have  resulted  in  a  new  class  of  compact,  ultrashort  pulse 
lasers  with  extremely  high  intensities  [1-3].  Intense  pulses  have  numerous  potential  applications 
in  areas  such  as  advanced  laser-driven  accelerators  [4-14],  harmonic  generators 
[15-19],  x-ray  lasers  [20,21],  other  short  wavelength  radiation  sources  [22],  and  “fast  ignitor” 
laser  fusion  [23-25].  Laser  technology  is  now  being  pushed  to  such  extremely  short  pulse 
lengths  that  the  pulse  is  only  a  few  optical  cycles  in  duration.  For  example,  at  a  wavelength  of 
0.8  Rm,  Baity  [2]  has  produced  pulses  of  4  TW  peak  power  with  durations  of  18  fs  (~  7 
wavelengths)  with  plans  to  extend  these  results  to  ~  10  fs  and  >  100  TW.  In  this  regime,  finite 
pulse  length  effects  may  play  an  important  role  in  the  laser  pulse  propagation  dynamics. 

The  propagation  dynamics  of  long  laser  pulses  (long  compared  to  the  wavelength)  is 
described  by  the  well-known  paraxial  wave  equation.  In  the  paraxial  wave  equation 
approximation,  lowest  order  diffraction  effects  associated  with  the  laser  beam  are  retained  but 
finite  pulse  length  and  higher  order  diffraction  effects  are  neglected.  The  solutions  to  the 
paraxial  wave  equation  in  vacuum  are  the  well-known  Laguerre-Gaussian  functions  that  describe 
the  dynamics  of  long  laser  beams  [26].  When  the  laser  pulse  length  becomes  sufficiently  short, 
i.e.,  less  than  ~  lO’s  of  wavelengths,  finite  pulse  length  effects  can  play  an  important  role  [27]. 
An  example  of  this  is  the  propagation  of  short  laser  pulses  in  a  guiding  channel.  Extended  laser 
pulse  propagation  in  a  plasma  channel  [28-40]  is  important  in  a  number  of  applications, 
including  high  gradient  accelerators  and  x-ray  lasers. 

In  this  paper  a  quasi  paraxial  approximation  (QPA)  is  introduced  which  is  an  extension  of 
the  well-known  paraxial  approximation  to  the  wave  equation  to  include  finite  pulse  length 
effects.  Employing  the  QPA,  a  pair  of  coupled  envelope-power  equations  is  derived  for  short 
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laser  pulses  propagating  in  vacuuni,  plasma  and  channels.  The  model  includes  atomic  electron 
and  relativistic  effects.  The  results  contained  in  this  paper  include:  i)  an  analytical  formulation 
of  short  laser  pulses  using  the  QPA,  ii)  a  derivation  of  a  pair  of  coupled  laser  envelope-power 
equations,  iii)  a  laser  envelope  modulation  which  eventually  damps  due  to  frequency  spread 
phase  mixing,  iv)  demonstration  of  significant  modification  of  nonlinear  processes  due  to  finite 
pulse  length  effects,  and  v)  analysis  of  short  pulse  propagation  dynamics  in  long  plasma 
channels.  These  results  represent  new  features  and  effects  associated  with  short  pulse  lasers. 
Short  laser  pulse  effects  have  been  eonsidered  numerically  or  in  one  dimension  by  others  [41- 
44],  however,  the  major  findings  in  this  paper  were  not  addressed.  For  example,  due  to  finite 
pulse  length  effects,  the  trailing  edge  of  an  unmatched  laser  pulse  propagating  in  a  plasma 
channel  is  found  to  undergo  an  envelope  oscillation,  while  the  leading  edge  is  damped.  For  a 
suffieiently  short  pulse  the  modulation  will  modify  the  main  body  of  the  pulse.  The  laser 
envelope  modulation  is  shown  to  be  due  to  the  dependence  of  the  pulse  group  velocity  on  the 
spot  size  through  the  pulse  length.  In  addition,  finite  pulse  length  effeets  are  shown  to 
significantly  increase  nonlinear  focusing  processes. 

The  organization  of  this  paper  is  as  follows.  The  general  wave  equation  for  the  electric 
field  of  a  finite  length  laser  pulse  propagating  in  vacuum,  plasmas  or  channels,  including 
nonlinear  (atomic  and  relativistic)  focusing  effects,  is  presented  and  discussed  in  Sec.  II.  In  Sec. 

III  the  propagation  dynamics  of  a  short  pulse  in  vacuum  or  uniform  plasma  is  derived.  In  Sec. 

IV  a  general  pair  of  coupled  envelope-power  equation  is  derived.  These  equations  are  used  to 
analyze  the  propagation  of  a  short  laser  pulse  in  a  preformed  plasma  channel  with  nonlinear 
focusing  effects.  Various  propagation  limits  are  discus.sed  and  the  laser  envelope  modulation  is 
analyzed.  The  laser  envelope  modulation  is  di.scussed  in  Sec.  V.  Numerical  illustrations  arc 
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presented  in  Sec.  VI  and  a  conclusion  is  given  in  Sec.  VII.  In  the  Appendix  the  range  of  validity 
of  the  QPA  is  obtained  in  two  limiting  cases. 


II.  Finite  Pulse  Length  Model 

The  propagation  medium  is  a  vacuum,  plasma  or  preformed  plasma  channel  consisting  of 
free  and  bound  electrons,  i.e.,  partially  stripped  plasma.  Nonlinear  processes  arising  from 
relativistic  and  atomic  polarization  effects  are  included.  The  wave  equation  includes  a  plasma 
current  consisting  of  free  electrons  and  a  polarization  current  arising  from  the  bound  atomic 
electrons  and  is  given  by  [31,45,46], 


where  is  the  transverse  Laplacian,  E(r,t)  is  the  electric  field,  Jp  is  the  plasma  current  density 
associated  with  the  free  electrons  and  P  is  the  polarization  field  associated  with  the  bound 
electrons.  The  atomic  polarization  field  consists  of  a  linear  and  nonlinear  part, 

P  =  (l/4;r)(77o  -  1  +  2j]Qri2l)^,  where  Tjo  is  the  linear  index,  r\2  is  the  nonlinear  refractive 
index  and  I  =  (c/47i)tio  (E  •  E)  is  the  time  averaged  laser  intensity.  In  the  present  model,  the 
origin  of  the  nonlinear  index  TI2  is  the  anharmonic  potential  in  which  the  bound  electrons 
oscillate.  Noting  that  dj p/dt  -  (4;r)“'a)p(r)E  and  setting  tjo  =  1,  the  wave  equation  becomes 


where  co,,  (r)=  (47i:qSip(r)/m)‘'"  is  the  plasma  frequency  and  np(r)  is  the  plasma  density.  In  Eq.  (2) 
the  coefficient  of  the  nonlinear  term  (3  =  Pp+  Pa  denotes  relativistic  free  electron  effects  as  well 
as  nonlinear  atomic  electron  effects  arising  from  T|2, 


me 


V  y 


(3a) 


o  ^  ® 


(3b) 


where  /m)*^^,  npo  is  the  density  on  axis  (r  =  0),  and  co  is  the  characteristic  laser 


frequency.  In  Eq.  (3),  Pp  results  in  a  critical  power  for  relativistic  focusing  [31,47-50]  while  the 
Pa  results  in  a  critical  power  for  atomic  electron  focusing  [31,51-53].  In  a  partially  stripped 
plasma,  atomic  effects  can  occur  on  a  time  scale  ~  lO  ’^  sec  and  can  dominate  free  electron 
effects  in  the  nonlinear  term  [32,45,46].  The  nonlinear  term,  due  to  relativistic  free  electrons,  is 
only  significant  when  the  laser  pulse  length  is  longer  than  a  plasma  period  [54].  The  critical 
powers  for  relativistic  focusing  in  plasma  and  nonlinear  focusing  in  a  gas  are,  respectively, 

Pp=lc{qlrY-{(iilO)p9,  (4a) 


and 


(4b) 


where  re  is  the  classical  electron  radius.  In  general,  when  the  laser  power  exceeds  either  of  these 
critical  powers,  focusing  occurs  [31,45,46].  The  total  nonlinear  focusing  power  consists  of 
contributions  from  both  Pp  and  Pa  and  is  given  by 

P,=PpPa/iPp+Pa)-  (4C) 
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The  preformed  plasma  density  channel  is  taken  to  be  a  function  of  radial  position  in  order 
to  provide  for  guiding  of  the  laser  pulse.  The  radial  dependence  of  the  plasma  frequency  is  given 
by 


(Op{r)  =  0}pQ 


1  + 


An 

npQ 


where  npo  +  An  is  the  density  at  the  edge 


xl/2 

(5) 


) 

of  the  plasma  channel  (r  =  Tc).  For  guiding,  the  plasma 


density  must  increase  as  a  function  of  r,  i.e.,  An  >  0. 

Laser  induced  plasma  waves,  wakefields  [55,56],  are  neglected.  This  can  be  justified  if 
the  laser  pulse  length  is  less  than  a  plasma  period  or  if  the  laser  intensity  is  sufficiently  low. 


a)  Quasi  Paraxial  Approximation  (QPA)  to  Wave  Equation 
The  laser  electric  field  is  of  the  form 

E  =  Eq  exp(/(fe  -  (ur))/  2  +  c.c.,  (6) 

where  Eo(r,t)  is  the  complex  amplitude,  k  is  the  wavenumber  and  (O  is  the  frequency. 
Substituting  Eqs.  (5)  and  (6)  into  Eq.  (2)  gives  the  wave  equation  for  Eo, 


Vi  +2/ 


0  £0  0 

0Z 


02^  0/2 


1  ^  -/i:2Z_  +  ^Eo-eS/2 


Eo  =  0,  (7) 


where  Kc  =  (a)po/c)(An/npo)*^^  is  the  focusing  parameter  associated  with  the  plasma  channel  and 
we  have  set  k  =  I -  (o],qIc^Y^ •  Changing  variables  from  (z,  t)  to  (z,^)  where  ^  =  z  - 


TipCt,  and  setting  k  =  TipOVc,  where  rip  =  (1  -  cojo  o^-axis  plasma  refractive 

index,  Eq.  (7)  becomes 
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V72  ^  ^ 

''1  +  2"!,,-^  + 2^^ 


+  (1-1^)-^ 
'  K 


+  :^-k2  — +  /;e„.e„/2 

9z  r/ 


Eo  =  0.  (8) 


The  second  term  on  the  left-hand  side  of  Eq.  (8)  represents  first  order  diffraction  effects,  the  third 
term  denotes  first  order  finite  pulse  effects,  the  fourth  and  fifth  terms  denote  higher  order  finite 
pulse  and  diffraction  effects  respectively  while  the  last  two  represent  guiding  and  nonlinear 
focusing,  respectively. 

In  the  absence  of  channel  guiding  (Kc  =  0)  and  nonlinear  focusing  (p  =  0)  we  can  obtain 
an  estimate  for  the  order  of  magnitude  of  the  various  terms  in  Eq.  (8).  The  second,  third,  fourth 
and  fifth  terms  in  Eq.  (8)  are  approximately  of  order 


2t] p{(i)lc)\dldz\  ~ 


.2  ’ 


T) 
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)4 
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respectively,  where  ro  and  are  the  spot  size  and  pulse  length. 


(9a) 

(9b) 

(9c) 

(9d) 

In  obtaining  the  estimates  in 


Eq.  (9)  we  used  |3/d^|  ~  IUq  and  |8/8z|  ~  l/Z/^,  where  Zr  =  TOipr^  IX  is  the  Rayleigh 

length  in  the  plasma  and  X  is  the  vacuum  wavelength.  The  relative  order  of  magnitudes  of  the 
first,  second,  third,  fourth  and  fifth  terms  in  Eq.  (8)  is 

i_  ^  .  A _ A 

Am]p  (q  ig  ■  AKijp  ■ 


1  :  I 


K1]p  t?o 


(10) 


The  first  two  terms  in  the  wave  equation  are  comparable  and  lead  to  the  paraxial  wave  equation 
approximation.  In  the  paraxial  approximation  finite  length  effects  and  higher  order  diffraction 
are  neglected,  i.e.,  the  third,  fourth  and  fifth  terms  in  Eq.  (8)  are  neglected.  The  paraxial  wave 

equation  assumes  that  XI Iq  «  1,  (1  -  liQ)XllQ  «  1,  and  X/Zr  «  1  and  is  given  by 

—  — ^0 

Solutions  of  the  paraxial  wave  equation  are  the  well-known  Laguerre-Gaussian  functions  [26]. 

When  the  Rayleigh  length  is  large  compared  to  the  pulse  length,  Zr  »  Iq,  the  higher 
order  diffraction  term  can  be  neglected  compared  to  the  finite  pulse  length  terms  resulting  in  the 
following  wave  equation 


f 


9  ft) 

V]_  +  2ir]  — 


^ _ i  c  3 

ft) 


az  ^  a^^ 


=  0. 


(12) 


Equation  (12)  contains  first  order  diffraction  and  finite  pulse  length  effects  and  reduces  to  the 
paraxial  equation,  Eq.  (11),  when  the  pulse  length  is  much  longer  than  the  wavelength,  Iq»X. 


b)  Finite  Pulse  Length  Wave  Equation 

Finite  pulse  length  effects  are  represented  by  the  terms  d  Idzd^  and  d  Id^  in  Eq. 
(12).  These  terms  under  certain  conditions  can  be  simplified  using  the  QPA,  allowing  for  the 
analytical  solution  for  the  field.  To  analyze  finite  pulse  length  effects  we  assume  that  the  general 
solution  of  Eq.  (12)  for  the  fundamental  transverse  Gaussian  complex  amplitude  is  given  by 

Eq  -  bexTp\(p  -  (\  +  i9)r^/ \ej_,  (13) 
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where  b,  (p,  0  and  are  real  functions  of  z  and  ^  =  z  -  rjpCt  and  is  a  unit  transverse  vector 

defining  the  polarization.  In  Eq.  (13),  b  is  the  amplitude,  (p  is  the  pha.sc,  0  is  related  to  the 
wavefront  curvature,  and  rj  is  the  spot  size  of  the  laser  pulse.  The  central  assumption  in  the  QPA  is 
that  the  main  contribution  from  the  finite  length  term  will  come  from  the  ^  dependence  in  the 
initial  amplitude.  Hence,  we  make  the  approximations  s  (3f/j(feo)/5'^)Eo 

where  bQ{^)  =  b(z  =  0,  ^).  In  the  Appendix 
the  QPA  will  be  shown  to  be  well  satisfied  for  a  broad  range  of  parameters.  Employing  the  QPA, 
Eq.  (12),  including  the  guiding  term  and  relativistic/atomic  electron  nonlinearities,  becomes 


3z 


En  =  0, 


where 


(14a) 


e(^)  = 


0)T]. 


d£n{bQ{^)) 

a^  ’ 


(14b) 


and 


/c\  c  de  2 

g(^)  =  - T7  -  e  • 

o)T]p  a^ 


(14c) 


Finite  pulse  length  effects  are  represented  by  the  functions  8(^)  and  g(^).  If  the  laser  pulse 


0  0 

amplitude  has  an  initial  Gaussian  longitudinal  profile,  ~  exp (-4^  /^o)’  ^hat 


e(^)  = 


8e  ^  4 


mp  £l  il' 


(15a) 


and 


8 


(15b) 


r  ^  f 

(N 

00 

1 

0 

1 

5(^)  =  2 

The  functions  £  and  g  have  magnitudes  much  less  than  unity  in  the  vicinity  of  the  laser  pulse,  i.e., 
1^1  <  f  0-  To  estimate  the  order  of  magnitude  for  e  and  g  we  take  X.  =  1  |im  and  f  q  =  ^5  )im  (50 


fs  pulse)  and  obtain  £(^  =  f  q)  =  8  x  10  ^  and  g(^  =  ^q)  =  -  6  x  10 

The  quasi-paraxial  wave  equation,  Eq.  (14a),  indicates  that  finite  pulse  length  effects  can 
be  important  when  fo  <  where  L  is  the  propagation  length.  This  implies  that  finite 

pulse  length  effects  are  significant  when  L»Z^,  as  would  be  the  case  for  a  guided  laser  pulse. 

The  laser  field,  which  is  described  by  Eq.  (13),  depends  on  the  functions  b,  (p,  0,  and  fs. 
To  obtain  equations  for  b,  9,  0,  and  rs,  we  substitute  Eq.  (13)  into  the  wave  equation,  Eq.  (14a), 
containing  finite  pulse  length  effects,  plasma  channel  guiding  and  nonlinear  focusing.  Equating 
like  powers  of  r  leads  to  coupled  equations  for  b,  (p,  0,  and  rs,[57] 


Q  +  i(\  +  ie)  — 


-£n(P/R^)  +  i(p 
2 


+P/(2R^)  -  a  =  0, 


az 


(16a) 

(16b) 


R 


where  Q(Z,^)  =  -(l  +  i6)/R^,  R  =  rs/ro  is  the  normalized  spot  size,  ro  is  the  initial  spot  size. 


Rm  =  rjro,  r^  =  {2vJKcf\  a(^)  =  (Z/j  /ro)^(l  -  Zr  =  TlpCOr^ /2c  =  ‘S  the 

Rayleigh  length,  Z  =  z/Zr,  P(Z,^)  =  Pl(Z,^)/Pc  is  the  laser  power  normalized  to  the  nonlinear 

focusing  power,  Pl(Z,^)  =  Ppe.-,k(^/^o)^  'S  the  laser  power,  Sq  is  the  peak  amplitude  at  Z  =  0, 

Ppeak  is  the  peak  laser  power  at  Z  =  0  and  Pc  is  the  nonlinear  focusing  power  in  Eq.  (4c)  which 
can  be  due  to  both  relativistic  and  atomic  electron  effects. 
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III.  Short  Pulse  Propagation  in  Vacuum  or  Uniform  Plasma 

a)  Fundamental  Transverse  Gaussian  Pulse  Solution 

In  this  section  we  obtain  and  discuss  the  finite  pulse  length  solution  to  Eq.  (14a)  by 
solving  Eqs.  (16a,b)  in  the  absence  of  guiding  (R,„  «x>)  and  nonlinear  focusing  effects  (Pc  — » 

««).  The  solution  of  Eq.  (16)  in  a  uniform  medium  with  refractive  index  rip  is 


b(Z,^)  =  bo(^)R-\z,^) 


\  +  £^iO 


l  +  e(^){Z  +  e{0) 


V/2  . 

exp 


l+£^((^) 


(17a) 


(p{Z,^)  =  -  tan-l(Z  +  £(^))  +  tan-J  (e(^))  -  — 4--  Z, 


(17b) 


e(z,^) 


(l  +  £(^)(Z  +  e(^))) 


/  .  ,^.\2 


R(Z,0  = 


1  +  (Z  +  e(^)Y 


1  +  £(^)(Z  +  £(0) 

and  the  wavefront  radius  of  curvature  is 

R,(Z,^)  =  r'^Zh  /e(Z4)  =  (Z^  /Z)(l  +  (Z  +  e(Z,Of). 


(17c) 


(17d) 


(17e) 


In  the  paraxial  limit,  £(^)  -4  0,  the  functions  b,  (p,  6,  R  and  Rc  in  Eqs.  (17)  reduce  to  the 
conventional  expressions  for  the  fundamental  transverse  Gaussian  beam  [26] 


b(Z)  =  bolRiZ) 

(18a) 

(p(Z)  =  -  tan“*  Z 

(18b) 

e(z)  =-z. 

(18c) 

r(Z)  =  (i+z2J  , 

(18d) 

R,{Z)  =  {Z^IZ){\  +  Z'^). 

(18c) 
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For  propagation  in  vacuum  we  set  T|p  =  1  in  Eqs.  ( 1 7)  and  ( 1 8). 


b)  Group  Velocity 


To  correctly  obtain  the  group  velocity  it  is  necessary  to  use  the  quasi-paraxial  wave 
equation  in  which  finite  pulse  length  effects  are  included.  The  paraxial  wave  equation  in  Eq. 
(11)  does  not  contain  finite  pulse  length  effects;  it  is  valid  for  infinitely  long  beams.  Since  Eq. 
(1 1)  is  independent  of  the  variable  ^  which  describes  the  pulse  shape,  finite  pulse  length  effects 
can  be  included  in  an  ad  hoc  manner  by  simply  multiplying  the  paraxial  solution  by  an  arbitrary 
pulse  envelope.  In  general,  the  group  velocity  is  the  velocity  with  which  the  peak  of  the  pulse 
envelope  travels,  i.e.,  the  velocity  for  which  db(Z,^)/d^  =  0.  In  the  paraxial  approximation  the 
group  velocity  is  found  to  be  =  cT]p  =  c{ckl(0)  where  we  have  taken  the  pulse  amplitude  to 


be  bQ{^)  =  bQ  exp(-4^^  /Iq)  and  Bq  is  the  peak  amplitude.  In  this  approximation  the  product  of 
the  group  and  phase  velocities  is  c  . 

In  the  quasi-paraxial  approximation,  which  includes  pulse  length  effects,  from  Eq.  (17a) 
the  pulse  envelope  can  be  written  as 

b{Z,0  =  ^o(^)t  -  /(I  +  zh\,  (19) 


where  bo  is  the  amplitude  in  the  paraxial  approximation  and  the  term  proportional  to  e  is  the  QPA 
correction  due  to  finite  pulse  length  effects.  Equation  (19)  is  correct  to  order  e  and  is  valid  for 
eZ  «  1 .  In  the  presence  of  finite  pulse  length  effects  the  group  velocity  is  found  to  be 


VgiO-^crii 


1-- 


^  A  ^ 


(l-cV/zj) 

(l^cV/Z^f 


^cr]p 


^  X  ^ 


717]  pro 


(20) 
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where  the  second  expression  is  valid  for  ct  «  Zr.  In  obtaining  Eq.  (20)  we  evaluated 


I 

dh/d^=0  in  Eq.  (19)  and  set  Z-ctlZj^  and  \dt'v^{t')-crjpt .  When  finite  pulse  length 

0 

effects  are  properly  included,  as  in  the  QPA,  to  lowest  order  in  e  the  group  velocity  is 
independent  of  the  pulse  length  but  depends  on  the  spot  size  as  given  in  Eq.  (20). 


IV.  Coupled  Envelope-Power  Equations 

In  this  section  we  analyze  the  propagation  of  a  finite  length  laser  pulse  in  a  preformed 
plasma  channel  including  nonlinear  focusing  effects.  In  general,  for  propagation  in  a  plasma 
channel  there  is  no  equilibrium  solution  for  a  finite  length  high  power  laser  pulse.  For  a 
sufficiently  low  power,  where  nonlinear  focusing  effects  can  be  neglected,  an  equilibrium 
solution  for  finite  length  pulses  exists.  However,  the  equilibrium  undergoes  an  envelope 
modulation.  The  normalized  envelope  equation  for  a  laser  pulse  in  a  plasma  channel,  to  first 
order  in  e(^),  is  obtained  from  Eqs.  (14a,b) 


where 


-  (1  -  P)R~^  +  R';;^R  +  eF{Z,^)  =  0, 


P  =  Po(^)exP 


£  JG(Z’,^)JZ’ 
0 


is  the  normalized  power, 

F{Z,^)  =  /?“^[3  +  {RIR,n)^  -  P+  {RdRldZ)^\)RldZ, 
C(Z,^)  =  y?"^[l  +  (/?//?,„)''  +  {RdRldZ)^\, 


(21a) 

(2Ib) 

(21c) 
(2  Id) 
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and  R  =  rs/io,  Z  =  z/Zr.  In  obtaining  Eqs.  (21)  higher  order  finite  length  effects  contained  in  a, 
which  are  of  order  e^  have  been  neglected.  The  coupled  envelope-power  equations  for  the 
guided  pulse  given  by  Eqs.  (21)  are  nonlinear  functions  of  spot  size  and  contain  finite  pulse 
length  effects.  Several  limiting  cases,  depending  on  the  laser  power  and  pulse  duration  will  be 
discussed.  In  the  following,  low  laser  power  refers  to  powers  much  less  than  the  nonlinear 
focusing  power,  Pl  «  Pc- 

a)  Low  power.  Long  pulse  (P  «  1,  e  =  0), 

In  the  low  power,  long  pulse  limit,  the  envelope  equation  reduces  to 

^  +  R-^R  -  R-^  =  0.  (22) 

3Z 

The  second  term  in  Eq.  (22)  denotes  plasma  channel  focusing  while  the  last  term  represents 
diffraction.  In  this  low  power,  long  pulse  limit,  the  beam  has  an  equilibrium  spot  size  Req  =  Rm- 

b)  Low  power.  Short  pulse  (P  «  1,  e  0), 

In  the  low  power,  short  pulse  limit  the  equilibrium  solution  to  Eqs.  (16)  is 


Req  Rm? 

(23a) 

beq(Z,^)  =  bo(^)exp|:-eZ/4j 

(23b) 

il 

P 

(23c) 

(23d) 

while  the  equilibrium  power  is  given  by 
=  Poexp[2e(p^.^{Z)], 
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(23e) 


and  the  pulse  group  velocity  is  given  by  the  expression  in  Eq.  (20).  The  effect  of  finite  pulse 
length  (e  ^  0)  is  to  make  the  peak  of  the  pulse  shift  backward  as  the  pulse  propagates. 


High  power.  Long  pulse  ( 1,  e  =  0) 

In  this  limit,  the  equilibrium  is  given  by 

(24a) 

(24b) 

(24c) 

where  Pgq  =  Po(^)  =  (Pmax/Pc)(^’o(^)^^o)  initial  normalized  laser  power  as  a  function  of  4- 

In  this  limit  the  laser  power  does  not  evolve  with  distance,  i.e.,  Peq  =  Po(^)  is  independent  of  Z. 
For  peak  laser  powers  less  than  the  critical  power,  Po(^)  <  1,  the  equilibrium  is  stable. 

In  this  limit  for  a  uniform  plasma  (no  channel  guiding,  Rm  =  ~),  a  uniform  laser  beam 
i.e.,  independent  of  will  have  a  longer  Rayleigh  length  longer  than  the  vacuum  Rayleigh 
length  because  of  nonlinear  focusing  effects.  The  modified  Rayleigh  length  is  given  by 

Zr-0-  Pe^r^'^Ziio,  (25) 

provided  Peq  <  1 .  For  Peq  >  I ,  the  beam  focuses  in  a  distance 

-ir‘'^Zso.  (26) 

In  Eqs.  (25)  and  (26),  Zro  is  the  Rayleigh  length  associated  with  low  power  lasers,  i.e.,  powers 
much  less  than  the  nonlinear  focusing  power  (Peq  «  I). 
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d)  High  power.  Short  pulse  (P<  1,  £  =  0) 


In  this  limit,  the  laser  pulse  does  not  have  an  equilibrium.  The  stability  characteristics  of 
laser  pulses  in  plasma  channels  can  be  obtained  by  perturbing  Eq.  (21)  about  the  spot  size  Req(4) 
=  Rm(l  -  Po(^))'^'‘,  where  R  =  Req  +  5R  and  i  5R|  «  Req.  The  equation  for  8R  is 

d^SRIdZ^  +  4R^ (1  - Po)6R  +  4ei?o^(l - PQ/2)d5R/dZ  =  -  Rq^6 P,  (27a) 


where 


f 

z 

SP  =  Po 

exp 

-e(GoZ  +  l5GdZ') 

-1 

V 

0 

J 

(27b) 


Go  =  2Ro  ^(1  -  Pq! 2)  and  5G  =  -  2Rq^Pq5  R.  Since  there  is  no  equilibrium  for  high  power, 

short  pulses  propagating  in  a  plasma  channel,  both  the  laser  power  and  spot  size  evolve  with 
propagation  distance.  The  spatial  evolution  of  the  power  and  spot  size,  among  other  factors,  is 
given  by  exp(-  2e  (1  -  Pq  /  2)Z  /  Rq  ) 


V.  Laser  Envelope  Modulation 

a)  Envelope  Oscillation 

In  the  low  power  (P  «  1),  short  pulse  limit  b)  of  the  preceding  section,  the  perturbed 
envelope  equation  is 


d^5R 

az' 


+  45R  +  4e 


ddR 

az 


^  c  d^6R  _ 

(0  a^az 


(28a) 


with  solution 


5R  -  SRq exp(-2fZ  -Z^  I z'^) cos(2Z) , 


(28b) 
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where  8Ro  =  8R(Z  =  0),  35R/3Z  =  0  at  Z  =  0  and  Zj-co(. ^J{2^[6c)  is  tlie  phase  mixing  length 
normalized  to  .  The  damping  term  is  obtained  from  an  extension  of  the  QPA  to  include  phase 
mixing.  In  obtaining  the  solution  in  Eq.  (28b)  it  was  assumed  that  the  envelope  period  kZ^  is 
small  compared  to  the  growth  length  Zr/(2I£|)  and  the  phase  mixing  length  (;i/6'^^)(^(/A)Zr  The 
modulation  amplitude  is  proportional  to 

exp[-8A^Z/(;r77/^„)  -  Z  VzJ] .  (29) 

The  modulation  at  the  front  of  the  pulse  (^>0)  is  always  damped  while  in  the  back  (<^<0)  it 
initially  grows  but  is  eventually  damped  due  to  frequency  spread  phase  mixing.  Overall 
damping  of  the  modulation  occurs  for  Z  >  2K\^\l{'iri pX) . 

b)  Laser  Modulation  Mechanism 

The  mechanism  for  the  envelope  modulation  can  be  understood  by  noting  the 
relationships  between  the  group  velocity,  spot  size  and  power  of  a  pulse  propagating  in  a  plasma 
channel.  The  group  velocity  of  a  pulse  in  a  plasma  channel  can  be  written  as  Vg  =  Vgo  +  6vg, 
where  the  mean  group  velocity  Vgo  is  given  by  Eq.  (20),  the  perturbed  group  velocity  is  5vg  = 
c(X/7tro)^8r/roand  5r  is  the  perturbed  spot  size.  To  lowest  order,  conservation  of  power  implies 
6b  =  -  bo6r/ro,  where  6b  and  bo  are  the  perturbed  and  unperturbed  laser  field  amplitudes, 
respectively.  Figure  l(a,b)  shows  the  amplitude  and  spot  size  of  a  finite  pulse  length  laser  in  a 
reference  frame  moving  with  the  mean  group  velocity  Vgo.  The  solid  curve  shows  the 
equilibrium  amplitude  and  spot  size  as  a  function  of  z  -  Vgot.  If  the  spot  size  is  uniformly 
increased  (6r  >  0)  along  the  pulse,  the  group  velocity  increases  by  the  amount  6vg.  The 
amplitude  in  front  of  the  unperturbed  pulse  increases  (6b  >  0)  while  the  amplitude  in  back  of  the 
unperturbed  pul.se  decreases  (6b  <  0).  Conservation  of  power  indicates  that  the  perturbed  spot 
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size  in  back  of  the  pulse  is  further  increased  since  6r  =  -  ro5b/bo  >  0  and  decreased  in  front  of  the 
pulse  since  6r  =  -  ro6b/bo  <  0.  If,  instead,  initially  the  spot  size  were  uniformly  decreased,  the 
spot  size  at  the  back  would  be  further  decreased  while  in  the  front  it  would  increase.  Hence,  the 
perturbed  spot  size  is  damped  in  front  of  the  pulse  and  is  unstable  in  the  back.  Substituting  ^  =  z 

-  Vgt  =  z  -  Vgot  -  6vgt  =  ^  -  6vgt  into  b  ~  boexp(-4^^/f  0 ).  the  rate  of  change  of  the  perturbed 
amplitude  is  ddbldt «  SbQ^odvg  /Iq.  Using  8r  =  -(r(/bo)5b,  we  find 

d5r/dt  =  -Sc^()iX,/7trQ)^Sr/£Q  which  agrees  with  the  growth  term  in  Eq.  (29).  Inherent  to  a 
finite  laser  pulse  is  a  frequency  spread  given  by  5(0  ~  d Hence  the  envelope  modulation 
frequency  =ldZ^  acquires  a  spread  5Qt~  ^  5(0/(0~Q.^{Xl27ii„) .  This  envelope  frequency 
spread  results  in  phase  mixing  of  the  modulation  in  a  distance  Zd  (normalized  to  Zr). 


VI.  Numerical  Illustrations 

Figure  2  is  a  plot  of  the  normalized  spot  size  R(Z,^)  in  Eq.  (17d)  as  a  function  of  ^  for 
various  values  of  Z  =  0,  1,  2,  3.  In  this  figure  the  laser  pulse  propagates  in  vacuum  and  the  pulse 
length  is  f  0  =  6  X.  In  the  absence  of  finite  length  effects  the  spot  size  would  be  independent  of 

Figure  2  indicates  that  the  tail  of  the  pulse  flares  out  more  than  the  front  of  the  pulse,  leading 
to  a  “trumpet”  pulse  shape. 

In  Figs.  3-6  the  laser  pulse  parameters  are  X  -  1  |im,  (0^(0  «  1  Iq  =20  p^m  (67  fs)  and 

the  peak  power  is  Ppeak  =  0.56  Pc.  The  total  nonlinear  focusing  power.  Pc,  is  given  by  Eq.  (4c) 
and  consists  of  contributions  from  free  and  atomic  electrons.  In  all  the  figures  there  is  an  initial 
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mismatch  in  (he  spot  size  compared  to  (he  equilibrium  spot  size,  i.e.,  Rq  =  R(0,0)  =  1  and  R,„  = 
1.15.  Figure  3(a)  shows  the  .spot  size  R(Z,^)  as  a  function  of  Z  =  z/Zr  and  IJX,  with  finite  pulse 
length  effects  (0  0)  included.  For  comparison  Fig.  3(b)  shows  the  same  plot  except  in  the 

absence  of  finite  length  effects  (e  =  0).  The  laser  envelope  modulation  is  clearly  seen  in  Fig.  3(a) 
where  the  spot  size  oscillations  at  the  front  of  the  pulse  (^  >  0)  are  damped  and  in  the  back  (^  < 
0)  grow.  Finite  pulse  effects  not  only  result  in  an  envelope  modulation  but  also  significantly 
enhance  nonlinear  focusing.  This  is  shown  in  Fig.  4  where  the  spot  size  with  finite  pulse  length 
effects  (solid  curve)  approaches  zero  at  =  -3  for  Z  =  15.  The  spot  size  without  finite  length 
effects  (dotted  curve)  shows  less  than  a  10%  decrease  at  ^  =  0  for  Z  =  15.  Figures  5(a)  and  (b) 
show  the  laser  pulse  amplitude  b(Z,^)  as  a  function  of  Z  and  ^X  with  and  without  finite  pulse 
length  effects,  respectively.  As  a  result  of  the  enhanced  nonlinear  focusing  due  to  the  finite  pules 
length.  Fig.  5(a)  shows  a  significant  increase  in  the  pulse  amplitude  at  Z  =  15  compared  to  Fig. 
5(b).  Figure  6  shows  the  pulse  power  as  a  function  of  ^X  and  radial  coordinate  r/ro  at  Z  =  15. 
Finite  length  effects  in  Fig.  6(a)  result  in  an  increase  in  the  peak  power  as  well  as  a  distortion  of 
the  pulse  compared  to  Fig.  6(b),  where  finite  length  effects  are  absent.  In  Fig.  6(a),  finite  length 
effects  reduce  the  pulse  propagation  velocity,  i.e.,  peak  of  the  pulse  occurs  at  negative  values  of 
In  Fig.  6(b),  nonlinear  focusing  effects  are  included  while  finite  pulse  length  effects  are 
neglected,  8  =  0.  For  e  =  0  the  pulse  velocity  is  c  and  nonlinear  focusing  is  substantially  reduced. 


VII.  Conclusions 

The  increasing  use  of  ultra  short  laser  pulses  in  many  applications  requires  that  the 
paraxial  wave  equation  be  extended  to  include  finite  pul.se  length  effects.  We  present  the  quasi 
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paraxial  approximation  (QPA)  to  the  wave  equation.  The  QPA  is  an  extension  of  the  usual 
paraxial  approximation  and  takes  finite  pulse  length  effects  into  account.  A  pair  of  coupled 
envelope-power  equations  is  derived  for  short  laser  pulses  propagating  in  vacuum,  plasmas  and 
preformed  plasma  channels.  The  model  includes  atomic  electron  and  relativistic  effects.  We 
find  that  finite  length  effects  can  significantly  modify  the  laser  field.  The  new  results  include:  i) 
an  analytical  formulation  of  short  laser  pulses,  ii)  a  derivation  of  a  pair  of  coupled  laser 
envelope-power  equations,  iii)  a  laser  envelope  modulation,  iv)  demonstration  of  significant 
modification  of  nonlinear  processes  by  finite  pulse  length  effects,  and  v)  analysis  of  short  pulse 
propagation  dynamics  in  long  plasma  channels. 
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Appendix:  Validity  of  the  Quasi  Paraxial  Approximation 

The  approximation,  which  leads  to  the  simplified  wave  equation  in  Eq.  (14a),  requires 
that 


£  En  » 


£  + 


T]„co 


Ef 


(Al) 


where  Eo  and  e  are  given  by  Eq.  (13)  and  (14b).  The  short  pulse  approximation  requires  that  Eq. 
(Al)  be  satisfied. 


a)  Approximation  in  Vacuum 

For  pulse  propagation  in  vacuum,  rip  =  1,  the  inequality  in  Eq.  (Al)  can  be  rewritten  and 
the  approximation  is  shown  to  be  valid  if 


Z 

1  +  i(Z  +  e) 


,  (rlrof 
1  +  i(Z  +  c) 


(A2) 


where  we  have  assumed  |  e|  «  1 .  The  finite  pulse  length  approximation  in  vacuum,  used  to 
replace  aEo/a<^  with  -  (o) / c)£(^)Eq  in  Eq.  (14a),  is  well  satisfied  everywhere  except  within  a 
small  function  of  a  wavelength  of  the  pulse’s  center. 


b)  Approximation  in  Guiding  Channel 

To  determine  the  validity  of  the  short  pulse  approximation  in  a  guiding  channel  we 
consider  the  low  power,  short  pulse  limit  of  Eq.  (21),  i.e.,  limit  b)  in  Sec.  IV.  The  field  in  Eq. 
(13)  for  the  equilibrium  solution  given  in  Eqs.  (23)  is 

E,^  =/>o(^)exp[-£Z/y?2 
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Substituting  Eeq  for  Eo  in  Eq.  (Al)  yields  the  condition  for  the  QPA  to  be  valid.  The 
approximation  is  valid  provided  |e|  »  (c/corjp) |9£/8^1z//?,?j,  which  for  a  pulse  having  a 
Gaussian  longitudinal  profile  is 

^\»-L±r;^z.  (A4) 

riplTt 

The  QPA  approximation  is  not  valid  for  long  propagation  distances. 
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Figure  Captions 


Fig.  1  Illustration  of  the  physical  mechanism  for  the  laser  envelope  modulation.  The  amplitude 
is  shown  in  (a)  and  the  spot  size  in  (b),  as  a  function  of  ^  =  z  -  Vgot,  where  Vgo  is  the 
mean  group  velocity.  The  solid  curves  correspond  to  the  equilibrium.  The  dashed  curves 
show  the  amplitude  and  spot  size  for  the  case  with  group  velocity  larger  than  Vgo  (1)  and 
the  case  with  group  velocity  smaller  than  Vgo  (2). 

Fig.  2  Plot  of  normalized  spot  size  R(Z,^)  as  a  function  of  ^  for  a  pulse  of  length  £q=6A, 

propagating  in  free  space.  The  four  curves  correspond  to  normalized  axial  points  Z  =  0 
(lowest  curve),  1,2,  and  3. 

Fig.  3  Surface  plots  of  spot  size  /?  as  a  function  of  <^  /  A  and  propagation  distance  Z  =  z/Z/j 
with  (a)  finite  pulse  length  effects  (8  ^  O)  and  (b)  finite  pulse  length  effects  neglected 
(e  =  O) .  The  parameters  are  A  =  Ijum,  £q=  20 fim,  Ppeak  =  0.56P^ . 

Fig. 4  Plot  of  spot  size  Pas  a  function  of  ^/X  after  a  propagation  distance  equal  to  15 
Rayleigh  lengths  (Z  =  15)  The  solid  (dotted)  curve  includes  (neglects)  finite  pulse 

length  effects.  Parameters  are  the  same  as  in  Fig.  3. 

Fig.  5  Surface  plots  of  laser  pulse  amplitude  b  as  a  function  of  ^  /  A  and  propagation  distance 

Z  =  z/Z/j  with  (a)  finite  pulse  length  effects  (fi  5^0)  and  (b)  finite  pulse  length  effects 
neglected  (e  =  O).  Parameters  are  the  same  as  in  Fig.  3. 

Fig.  6  Surface  plots  of  laser  pulse  power  P  as  a  function  of  and  radial  coordinate  r/ro  after  a 
propagation  distance  equal  to  15  Rayleigh  lengths  (Z  =  15).  In  (a)  finite  pules  length 
effects  are  included  and  show  enhanced  focusing  and  decreased  propagation  velocity,  i.e.. 


peak  of  the  pulse  oceurs  for  negative  values  of  4-  In  (b)  finite  pulse  length  effeets  are 
neglected,  i.e.,  e  =  0  while  nonlinear  focusing  effects  arc  included.  Parameters  arc  the 
same  as  in  Fig.  3. 
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